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Equilibrium magnetization distribution of the vortex state magnetic nanoparticle is affected by the influence 
of the spin-transfer torque: an off-center out-of-plane vortex appears in the case of the disk shape particle and 
pure planar vortex in the case of asymmetric ring shape particle. The spin current causes extra out-of-plane 
magnetization structures identical to well known dip structures for the moving vortex. The shape of the dip 
structure depends on the current strength and value of the off-center displacement and it does not depend on 
the vortex polarity. The critical current depends on the nanodot thickness. 

PACS numbers: 75.10.Hk, 75.40.Mg, 05.45.-a, 85.75.-d 



The control of magnetic nonlinear structures using an 
electrical current is of special interest for applications 
in spintronicsP^ The spin-transfer torque acts on nonho- 
mogeneities in magnetization distributions, in particular, 
on magnetic vorticesP^ In particular, it w as p redicted 
theoretically and confirmed experimentally^' that the 
vortex core magnetization (so-called vortex polarity) can 
be switched on a picosecond time scale. This discov- 
ery demonstrates the potential of realizing all-electrically 
controlled magnetic memory devices, changing the direc- 
tion of the modern spintronicsP Recently Shibata et alW 
used the in-plane spin current to demonstrate the effect 
of the spin-transfer torque on the vortex state magnetic 
nanodisk. The spin current excites the spiral motion of 
the vortex which finally relaxes to some shifted position. 
Such a picture is a result of a Thiele-like vortex dynam- 
ics, where the vortex does not change the shape during 
its motion and it is valid for the relatively small vortex 
shifts. 

In this Letter we predict the dip formation nearby the 
vortex under the influence of the current. We consider 
two geometries, nanodisk and asymmetric nanoring. The 
spin current in a latter case provides the pure planar vor- 
tex centered on inner hole of the ring. Using micromag- 
netic simulations^ we discovered that such a shifted im- 
mobile vortex gets an extra out-of-plane magnetization 
which corresponds to the well-known "dip" structure of 
the moving vortexP Our study shows that the dip de- 
velopment is specified by the current direction and its 
intensity, and it does not depend on the vortex polarity. 
We also show that the critical current strongly depends 
on the nanodot thickness. 

Our study is based on a simulated magnetization dy- 
namics in the framework of the modified Landau-Lifshitz 
equation with the adiabatic spin torque term^ 

rn = —jm x H e g + am x rh — (u ■ V) m. (1) 
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FIG. 1. Profiles of the immobile equilibrium vortices for dif- 
ferent values of the applied current. The profiles were taken 
along the diameter line perpendicular to the current direction 
- dot-dashed line on the insets a) and b). The insets demon- 
strates the 2D distribution of the out-of-plane magnetization 
for two cases: zero current - inset a) (the corresponding pro- 
file is shown with thin solid line) , and current close to critical 
one - inset b) (the corresponding profile is shown with thick 
solid line). Dashed isolines corresponds to the case m z < 
and solid isolines - to m z > 0. All presented data were ob- 
tained from the micromagnetic simulations of the permalloy 
nanodisk with radius L — lOOnm and thickness, h — 20nm. 



Here m is normalized magnetization vector: m = 
M/Mg with Ms being the saturation magnetization, 
7 > is giromagnetic ratio, H e ff is the effective micro- 
magnetic field, a is the Gilbert damping constant. The 
velocity u is directed along the direction of electrons flow, 
with an amplitude u oc j, where j is the current density, 
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FIG. 2. Parameters of the equilibrium vortex state depend- 
ing on the applied current. Plot a) shows the vortex (solid 
line) and dip (dashed line) displacement from the disk center 
depending on the velocity u. The displacements are measured 
in units of L. Plot b) shows dependence of the dip depth rrid 
on u: points represents the simulation data and curve cor- 
responds to Eq. |2|. The vertical dashed line denotes the 
critical current u c when the vortex-antivortex pair is born. 
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FIG. 3. In-plane vortex profiles for different current values 
and vortex positions. Dashed and solid lines corresponds to 
the opposite current directions. The inset a) shows profiles 
of the initially in-plane vortex pinned by the hole displaced 
by the value L/2 from the disk center. The inset b) shows 
the similar profiles for the centered vortex. Simulations were 
performed for the permalloy disk with L — 100 nm and h = 10 
nm. 



which is supposed to be spatially uniform and constant. 
The permalloy^ disk with radius L = lOOnm and thick- 
ness h — 20nm was chosen for the numerical experiment. 
The ground state of such a disk is the out-of-plane vor- 
tex, situated at the disk origin [see Fig. [lp,)]. Then the 
spin current velocity was switched on adiabatically: step- 
by-step we increased the value of velocity u by a small 
value. On each step the full relaxation was achieved and 
the set of parameters was determined: position of the 
vortex center s, position of the dip minimum Sd and also 
the dip depth md (see Fig. [I]). As one can see Fig. 2j the 
vortex and dip displacements linearly depend on the ap- 
plied current. When the vortex shifts the dip amplitude 
increases and dependence rrid(u) is essentially nonlinear, 
[see Fig.[2j inset b)]: 
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when u <C u c , 
1 when u < u c . 
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When the velocity u achieves some critical value u c 
the dip structure become unstable and depth abruptly 
achieves its minimum value = — 1, hence the vortex- 
antivortex pair is born. Then the vortex polarity switch- 
ing occurs accordingly to the well studied scenarioP 
It should be noted that the obtained critical velocity 
u c = 355m/s is very close to the critical velocity for 
moving vortexP^ 

In order to study the role of the vortex out-of-plane 
structure in the dip formation process we performed the 
second kind of numerical experiment. With this end in 
view we considered the shifted vortex without out-of- 
plane component, the pure immobile planar vortex. The 
vortex was pinned on the small hole placed in the half 
of the disk radius. The hole radius r/, was chosen to be 
in the range € (a c ;r c ~ £), where a c is the critical 
radius of the transition between out-of-plane and planar 
vortices, 13 r c is vortex core radius and £ is the exchange 



length 11 . Thus the hole radius r/, is big enough to make 
the vortex planar and it is small enough to prevent signif- 
icant influence on the out-of-core magnetization struc- 
tures. For thin permalloy films a c — 0.37£ « 2 nnf^ and 
^ ~ 5.3 nm 11 . That is why we chose = 4 nm. For sim- 
ulation was chosen disk with L = 100 nm and h = 10 nm. 
The thickness decrease is needed for increasing the pin- 
ning effect due to decreasing the volume magnetostatic 
charges. 

Let us consider the planar vortex which is pinned at 
the inner hole of the ring. In the case of the symmetric 
hole, the planar vortex forms a ground state of the disk 
[see Fig. [3}})]. In the case of the asymmetric hole, the 
planar vortex is a equilibrium state, which corresponds 
to the local energy minimum in the current absence. One 
can see from the Fig. [3] the vortices does not have core 
or any other out-of-plane components for the case of 
zero current. However under the influence of the current 
there appears a dip-like out-of-plane structure. This 
dip structure has the following main properties: (i) The 
dip structure of the centered vortex is symmetrical and 
its sign is changed when the sign of u is changed, see 
Fig. ^ijp) . (ii) The skewness appears when the vortex is 
shifted and its value rapidly increases when the displace- 
ment s increases, (iii) The sign of the dip structure of the 
shifted vortex is determined by the sign of the product 
su and it has no direct relation to the vortex polarity, 
(iv) The critical current u c takes maximal values for the 
centered vortex and decreases as the vortex displacement 
s increases. The detailed analytical description of the dip 
development, based on the magnon mode analysis, is un- 
der consideration. 

When the effective current velocity u achieves its crit- 
ical value u c the vortex-antivortex pair is born on the 
edge of the inner hole. The antivortex falls into the hole 
(it annihilates with the pinned vortex) , so finally the new 
vortex appears out of the hole. Polarity of the new born 
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FIG. 4. Dependencies of the dip depth on the applied current 
for different thicknesses. Inset a) - in-plane vortex is pinned 
on the hole with s — 0.5, b) - free out-of-plane vortex. Solid 
black lines demonstrates the dependence The critical 

velocity u c = 330 ± 37 m/s which was determined in Ref. Q2] 
is shown in the subfigure b) as vertical strip. 

vortex p is equal to the polarity of the dip, in the other 
words p = — sign(us). In the case of centered vortex 
(s = 0) the polarity p is determined in a random way. 

It is interesting to note that the critical current value 
depends on disk thickness and this dependence is differ- 
ent for free (out-of-plane) and pinned (planar) vortices, 
see Fig. |4| In the case of free vortex the critical velocity 
u c weakly dependents on thickness while in the case of 
pinned vortex this dependence is drastic. 

In conclusion, we found that immobile planar vortex in 
a nanoring forms a dip structure under the action of adi- 
abatic spin current under the threshold value. Above its 
value the vortex-antivortex pair can be nucleated, which 
can cause the vortex switching process. 
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